Intrinsic coupling of orbital excitations to spin fluctuations in Mott insulators 
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We show how the general and basic asymmetry between two fundamental degrees of freedom 
present in strongly correlated oxides, spin and orbital, has very profound repercussions on the 
elementary spin and orbital excitations. Whereas the magnons remain largely unaffected, orbitons 
become inherently coupled with spin fluctuations in spin-orbital models with antiferromagnetic and 
ferroorbital ordered ground states. The composite orbiton-magnon modes that emerge fractionalize 
again in one dimension, giving rise to spin-orbital separation in the peculiar regime where spinons 
are faster than orbitons. 
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In transition metal oxides, different 3d orbitals near 
the Fermi level can have similar energy and thereby con- 
tribute to the low-energy physics. In presence of strong 
correlations charge fluctuations become suppressed and a 
Mott insulator is realized when there is a commensurate 
number of electrons per unit cell. The effective Hamilto- 
nian that emerges, often referred to as Kugel-Khomskii 
(KK) Hamiltonian, can be expressed in terms of spin and 
orbital operators [if. From a formal viewpoint the spin 
and orbital operators are very similar because they form 
identical algebras, but the way in which these operators 
enter into realistic KK Hamiltonians is very different. 

While the spin wavefunction is in essence rotationally 
invariant and the coupling between spin operators there- 
fore SU(2) symmetric, the symmetry in orbital space 
is much lower due to the ubiquitous crystal field act- 
ing on the orbital wavefunctions in crystals. This re- 
duced symmetry causes the orbital-orbital interaction to 
be very anisotropic in space and often inherently frus- 
trated, which causes exotic effects such as macroscopic 
degeneracy of the groundstate or the emergence of non- 
Abelian topological excitations in the case of the compass 
or Kitaev models, respectively |2j , which may be relevant 
for quantum computation [3, We will show that this 
asymmetry between spin and orbital degrees of freedom 
has fundamental repercussions on the coupling of the 
elementary spin and orbital excitations - the magnons 
and orbitons, respectively: whereas the magnons remain 
largely unaffected, orbitons become inherently coupled 
with spin fluctuations in spin-orbital models with anti- 
ferromagnetic (AF) and ferroorbital (FO) ordered ground 
states. This is relevant in the experimental context 
as substantial progress is being made in measuring or- 
bitons 0,0] and their dispersion, inparticular in resonant 
inelastic x-ray scattering (RIXS) [JS], where we predict 
the coupling of the orbiton to magnetic fluctuations to 
be clearly discernible. 

In the standard approach the complex problem of in- 
tertwined spin-orbital excitations is solved using a mean- 
field decoupling, i.e. considering magnons in a fixed or- 



bital background or orbitons in a fixed spin background. 
While such an approach can work well to obtain the 
correct spin and orbital orderings consistent with the 
Goodenough-Kanamori rules [l| and in some cases with 
ferromagnetic (FM) order we show here that it fails 
to describe orbital excitations even qualitatively correctly 
for a number of spin-orbital models. We avoid this decou- 
pling by mapping the coupled orbiton-magnon dynamics 
onto the well-controlled problem of a hole propagating 
in a magnetic background: the extensively studied sin- 
gle hole t-J model. In particular, we find that in one 
dimension (ID) the orbital excitation fractionalizes into 
freely propagating spinon and orbiton, giving rise to spin- 
orbital separation in the peculiar regime where spinons 
are faster than orbitons. 

Model and problem statement. — The generic form of 
the KK Hamiltonian [l[ in the Mott-insulating limit is 



(1) 



d,i) 



where i and j are lattice sites and on each bond the 
spin-spin interaction is Hfj whereas the orbital-orbital 
one is Hj^, To capture the generic differences between 
orbitons and magnons, it is enough to break the SU(2) ro- 
tation symmetry in orbital space, which is here achieved 
by the large local crystal field breaking the degeneracy 
between the orbitals, expressed as Hf = §fTf. We 
keep for simplicity the rotational symmetry in the in- 
teractions so that for the spins Hf- = S{ ■ Sj + \ and 
for the orbitals Hf^ = T,- ■ Tj + |, where S (T) arc the 
spin (orbital) operators that fulfill the SU (2) algebra for 
S = 1/2 (T = 1/2) spins (pscudospins) . The constant 
J > gives the energy scale of the spin-orbital superex- 
changc and the symmetry breaking field for the orbitals 
is E z . Note that for the case of E z — (not considered 
here), this model has an SU(A) symmetry even higher 
than the combined SU(2) x SU(2) symmetries, which re- 
sults in the ground state given by the Bethe Ansatz and 
composite spin-orbital gapless excitations in addition to 
the separate spin and orbital ones (see, e.g., Ref. fiol]). 
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This model describes the low-energy physics, determined 
by singly occupied sites, of a two-orbital Hubbard model 
in the limit of a large onsite Coulomb repulsion U and 
vanishing Hund's exchange Jh, cf. Eqs. (1-5) and (29) 
in Ref. [9J. 

Here we are interested in the orbital excitations of the 
model Eq. ([lj when E z 3> J i.e., the orbital splittings are 
larger than magnetic coupling energy. This is a realistic 
regime for many strongly correlated compounds such as 
ID or two-dimensional (2D) cuprates As one degree 
of freedom is completely polarized, the ground state is 
easily found and given by all electrons occupying a single 
orbital in an AF state. But while decoupling spin and or- 
bital degrees of freedom works here for the ground state, 
it is not at all appropriate for orbital excitations - this 
is the problem investigated below. 

Decoupling of spin and orbital sector. — The ground 
state = \ips) ® IV'o) of Eq. (QJ is described by the 
ground state \ips) = |AF) of an AF Heisenberg system 
formed by spins in the lower-energy orbital, i.e., an FO 
ordered state \tpo) = |FO). The orbital excitations are 
reached by flipping an orbital, i.e., by promoting an elec- 
tron at site j from the occupied lower orbital to the empty 
higher band at the same site, expressed by the orbital 
raising operator . The momentum-dependent orbital 
excitation is given by T+ = £\ e^T+ [T fc _ 
and the spectral function describing its dynamics is 
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First we discuss the orbiton spectral function by taking 
the orbital excitation to be independent of the magnetic 
excitation. One can then rewrite the orbital operators 
by use of Holstein-Primakoff bosons (see, e.g., Ref. [HJ]), 
keep only quadratic terms in the expansion (orbital- wave 
theory) and, noticing that the ground state \i/jo) does not 
contain bosons, obtain 



0(k,u) = 5[lj — LJow(k)], 



(3) 



with a mean- field orbital- wave dispersion uow{k) = 
E z — izJow(l — Ik)- Here, z is the coordination num- 
ber, 7^ is the lattice structure factor, and the effective 
orbital exchange constant Jow = 2J(V>s|Si • Sj + zlV's)- 
The orbital excitation on the mean-field level is thus a 
quasiparticle with a cosine-like dispersion with period 
2tt: for example in ID we obtain an effective reduced 
Jow — —0.4 J <C J, cf. the thick line in Fig. [2j An anal- 
ogous procedure for magnons in FM planes with alternat- 
ing orbitals (AO) has been applied to LaMn03 or KCUF3, 
and similarly yields magnons with a reduced bandwidth, 
but without any other trace of the AO order, in agree- 
ment with experiment 



13|. 



We will see, however, that 
for orbitons this framework of mean-field decoupling is 
greatly oversimplified. 

Mapping onto an effective t-J model. — The orbital- 
wave approximation puts all the impact of the AF order 
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FIG. 1: (color online) Two superexchange processes moving 
an electron in an excited orbital (indicated by oval) from site 
j to its neighbor j + 1: (a) and (b) describe orbiton motion 
when spins along the bond are antiparallel or parallel, re- 
spectively, see text. The states in the grey middle panels are 
not part of the low-energy Hilbert space corresponding to fl} . 
These virtual excitations within the full two-orbital Hubbard 
model illustrate the origin of those superexchange interactions 
of Hamiltonian (ft) that propagate the orbiton; note that the 
spin of the excited electron is conserved. 



into Jow: which is a mean-field average of the sum of 
two distinct superexchange processes in which the orbital 
excitation may propagate through the lattice. The first 
one, which corresponds to (T^T~ +1 +h.c.)(S^ Sf +1 +h.c.) 
processes in Hamiltonian ([T]), allows for orbiton prop- 
agation when the spins on the bond are antiparallel, 
see Fig. 0Ja). The second one, which corresponds to 
(TfTf +1 + h.c.)(S*S* +1 + 1/4) processes in Hamiltonian 
(PJ , allows for orbiton propagation when the spins on the 
bond are parallel, see Fig.QJb). Crucially (see next para- 
graph), this figure illustrates that in both cases the spin 
of the electron in the upper orbital 2 is conserved dur- 
ing the orbiton propagation. This is because spin-orbital 
Hamiltonian ([!} is, as mentioned above, a low energy 
limit of the two-orbital Hubbard model with the Hund's 
exchange Jjj — and the spins of individual electrons in 
the superexchange process cannot be flipped (see middle 
panels of Fig. [T]). In more realistic spin-orbital models 
the Hund's exchange is finite [lj], but the processes which 
would not conserve the electron's spin in the excited or- 
bital are small (oc Jh/U) and thus could be neglected. 

To be explicit, we now focus on the ID case and em- 
ploy a Jordan- Wigncr transformation [14| . 2D and three- 
dimensional (3D) cases are discussed afterwards. We 
thus introduce 5+ = fte" Q , Sf = e~™ Q f3], S*| = 
rij/3, where Q = J^ti n W an< ^ Pj create spinons 
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while T+ = c-™®a] 
where Q = £)f=i riu 



and 
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U 3 S ! j ~ "J" 2' 

creates a pseudospinon. 
Since the spin of the propagating electron in the up- 
per orbital 2 is conserved (see above), one may cal- 
culate the spectral function 0(k,cj) for, e.g., spin-up 
in the upper orbital. We are then allowed to replace 
T l -> Ej c lk3 T+{\ + SI) = c-"«a](l - n j0 ) in Eq. © 
and thus terms in the Hamiltonian that create or annihi- 
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late both a spinon and a pseudospinon at the same site 
lead to a vanishing contribution to the spectral function. 
Phase factors cancel in one dimension and we obtain 

0(ft,M) = ilim3(#» fc * , a+hft), (4) 

7r 7;->o u+E^—n—E z ~ir] 

with the effective fermionic Hamiltonian 



j+h.c.)--(ni0+n j p)+n i pri j j3 , (5) 
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and an implicit constraint Vj + a^ay < 1. Here |V>) 
is a tensor product of the magnetic ground state \ips) = 
|AF) expressed in terms of spinons and a vacuum state 
for pseudospinons [recall that we consider here a single 
orbiton only; this also allowed us to skip quartic terms 
in pseudospinons in Eq. ([15]) ]. 

At this point, we observe that the resulting effective 
Hamiltonian is in fact a Hamiltonian for the t-J model 
written in terms of the Jordan- Wigner fermions with the 
above constraint [la . By introducing the electron oper- 
= ajfije™® acting in the restricted 
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Hilbert space without double occupancies we obtain 



0(k, oj) = — lim 

7T >)^0 
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with the t-J Hamiltonian 



-Er-n-E z -iT] 



Pkt\4>), ( 6 ) 



H- 



whcr 



(p\<jPi° + h - c -) + J H^ Sl ■ S i + 7^ n J")' ( 7 ) 
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a and the hopping parameter t is 
defined as t — J/2 [16(. The ground state \ip) is now the 
tensor product of a vacuum state for holes and the \ips) 
state. We have thus mapped the single orbiton in the FO 
and AF chain, with dynamics governed by Hamiltonian 
([I]), onto a single hole doped into the undoped AF chain 
with its dynamics governed by Hamiltonian (|17l) . 

Numerical results for the t-J model. — To flesh out 
the resulting coupling between orbitons and spin fluc- 
tuations we use Lanczos exact diagonalization to evalu- 
ate Eq. (fTB|) on a finite chain (28 sites). The spectral 
function is shown in Fig. [2j the spectrum differs qualita- 
tively from the orbital-wave result shown as a thick line 
in Fig. [21 It now consists of multiple peaks (expected to 
merge into incoherent spectrum in the thermodynamic 
limit) instead of one single excitation. There is a dom- 
inant feature at the lower edge of the spectrum, but its 
periodicity is ir reflecting the doubled unit cell of the AF 
order. 

The t-J model with J > t [l6[ is not easily accessi- 
ble in the Hubbard type models, as it would formally 



FIG. 2: (color online) Spectral function 0{k,ui) of orbital 
excitation obtained via the mapping onto the t-J model, Eq. 
(|16[) . evaluated using Lanczos exact diagonalization on a 28 
site chain. A broadening r\ = 0.03J and E z = 10J. The thick 
line shows orbital excitation in a mean-field (orbital-wave) 
approach, Eq. 
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FIG. 3: (color online) Schematic representation of the orbital 
motion and the induced spin fluctuations giving rise to spin- 
orbital separation in ID. The first hop of the excited state 
(a —¥ b) creates a spinon (wavy line) that moves via spin 
exchange oc J. Next hop (b — > c) does not produce any extra 
spinons: an 'orbiton' freely propagating as a 'holon' with an 
effective hopping t ~ J/2 is created. 



correspond to small onsite interaction U, where the t— 
J model is no longer valid. In this regime, the spinon 
moves faster than the holon, and the entire lower edge 
of the spectrum is thus given by 'holon' states [It}- If 
the orbiton takes the place of the holon as argued here, 
this exotic behavior should be observable in RIXS exper- 
iments and one would expect a dominant excitation with 
orbiton-character with dispersion uj rj E z — 2tsin|fc| at 
the bottom of the spectrum. The latter would be ex- 
pected to extend up to cj w E z + \/J 2 + At 2 — At J cos k 
and to contain an intermediate feature still well-visible 
within the continuum with the dispersion of the purely 
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orbiton-character scaling as uj sa E z + 2tsin \k\. 

Figure [3] illustrates how electron exchange processes 
can let an orbital excitation propagate through the sys- 
tem after creating a spinon in the first step. The spinon 
itself moves via spin flips oc J > t, faster than the or- 
biton, and the two get well separated. The orbital-wave 
picture, on the other hand, would require the orbital ex- 
citation to move without creating the spinon in the first 
step. As can be inferred from Fig. 02 this is only possible 
for imperfect Neel AF spin order so that the averages of 
processes shown in Fig. [Ua) and (b) are finite. 

2D and 3D. — Remarkably, the standard OW picture 
becomes even worse in higher dimensions: in 2D (3D) cu- 
bic lattices, the mean-field orbital coupling Jow almost 
vanishes due to {ip s \Si ■ $j + \\^s) - -0.08 (-0.05); the 
orbital dynamics is then entirely governed by coupling 
to the spin fluctuations [l8|. The t-J model description 
of a single orbiton Eqs. (| 16411 T[) is general and valid for 
any dimension: the mapping rests entirely on the fact 
that spin on the excited orbital is conserved. In particu- 
lar, the Jordan- Wigner fermionization applied to the 2D 
case gives the t-J model expressed in terms of Jordan- 
Wigncr fcrmions II 511 and consequently the 2D version of 
Eqs. (jl6H17p . see ]ig| for details. 

A very good approximate solution for the higher- 
dimensional t-J model with J > t can be obtained 
by perturbation theory, because J > t corresponds to 
weak coupling, see Eq. (7) of Ref. (20|. The solu- 
tion shows that in 2D or 3D a single orbiton in the 
undoped AF system described by the t-J model can- 
not fractionalize due the magnetic string effect. Still, 
the orbiton dressed with spin fluctuations is mobile on 
a renormalized scale [23]. For example the 2D case di- 
rectly corresponds to the result in Ref. [2l| which gives 
w~£ z -ii +£2 (cos k x + cos k y ) 2 + xz(cos2k x + cos 2k y ) 
for the orbiton dispersion (xi are positive parameters 
cx t 2 j J when J > t. cf. Ref. [2l| for exact values). 
We thus expect to observe, e.g. in high resolution RIXS 
experiment on 2D cuprates, a spectrum similar to the 
one-particle Greens function of the higher-dimensional 
t-J model, with a low-energy quasi- par ticle orbiton peak 



cx bj and additional incoherent part [2l[ . Crucially, while 



the parameters Xi depend on J/t and would vary in more 
realistic spin-orbital models [16| , the general shape of the 
orbiton dispersion w - the minimum at (7r/2,7r/2), sad- 
dle point at (tt, 0), maxima at (0, 0) and (-7T, tt) - is robust 
reflecting the fact that the coherent motion of orbiton is 
possible only within a given spin-sublatticc. 

No analogue in spin sector. — Let us now revisit the 
analogue where the spin is polarized (FM) and the orbital 
sector shows AO order due to, e.g., a large J aim- Teller 
effect. An example would be the FM and AO planes in 
LaMnOa or KCUF3 Many studies have shown that 
the orbital degrees of freedom merely renormalize the 
spin excitation in this case and do not change the peri- 
odicity of magnons (see, e.g., Ref. [22|). This qualitative 



difference is caused by the fact that in realistic cases the 
Jahn- Teller stabilized AO order is much more classical 
and robust, thus suppressing the creation of orbital ex- 
citations by the excited spin and giving larger weight to 
the pure spin excitation. In other words, spin waves are 
typically below the orbital gap and well protected by the 
underlying SU(2) symmetry and Goldstone theorem. 

In conclusion, we have shown that orbitons in realistic 
spin-orbital models with AF-spin and ferroorbital ground 
state are so strongly coupled to the spin excitations that 
the usual mean-field decoupling of two sectors breaks 
down. In fact, we have presented an exact mapping of 
the problem onto an effective t-J model and have shown 
that the orbiton in such models behaves like a single hole 
in undoped AF Mott insulator. However, since the typi- 
cal supcrcxchangc parameters for KK spin-orbital models 
lead to J > t, the study of orbiton problem provides ac- 
cess to a regime of the t-J model which has been thought 
to be 'unphysical' in terms of single-band Hubbard mod- 
els. Finally, in ID signatures of spin-orbital separation 
are expected again in the peculiar regime where spinons 
are faster than orbitons. 
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SUPPLEMENTARY MATERIAL 

In what follows we show that the t-J model descrip- 
tion of the orbiton problem [Eqs. (6-7) in the main text] 
is not only valid in the ID case but also in higher di- 
mensions. To be explicit we concentrate now on the 2D 
case (from which the 3D case follows in a straightforward 
way) and introduce the Jordan- Wigner fermions a and /3 
for pseudospins and spins: 



St = pjj* Q i , 
S7 = e-™ Q ipl 



Si 



where Qj = Xw = i n W anc ^ 01 create spinons while 



TP 



17 = Ot^ Qi 

J J 

Ta — 71 -in. ~~ 



(8) 
(9) 
(10) 
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where Qj = X/i=i n '"- Note that this is the same trans- 
formation as in the main text but, to keep track of the 
phase factors, we explicitly wrote the site index j of the 
phase factors Qj and Qj. 

Next, similarly as in ID, the spin of the propagating 
electron in the upper orbital 2 is conserved, and one may 
calculate the spectral function 0(k,uj) for, e.g., spin-up 
in the upper orbital. We are then allowed to replace 



T l -> Hi e^T+(i + S Z A = e-^aUl 



n 3P) 



Eq. (2) 



in the main text and thus terms in the Hamiltonian that 
create or annihilate both a spinon and a pseudospinon at 
the same site also lead to a vanishing contribution to the 
spectral function. 

In 2D one has to take care of the phase factors. How- 
ever, the crucial observation is that for the case of the sin- 
gle orbiton the phase factors Qj associated with a pseu- 
dospinon either do not contribute at all [Eq. (|14p below] 



or cancel for the nearest neighbor bonds [Eq. (fl~5|) below] 
- similarly to ID. Thus, these are only the spin phase 
factors Qj which are in the end present in the spectral 
function and Hamiltonian written in terms of the spinons 
and pseudospinons: 



0(k,u) = - lim3( 

7T l)->0 
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with the effective fcrmionic Hamiltonian 
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As stated above (and mentioned for the ID case in the 
main text), here we have an implicit constraint Vj 0\Pj + 
ctjCtj < 1 while is a tensor product of the magnetic 
ground state \ips) = |AF) expressed in terms of spinons 
and a vacuum state for pseudospinons. 

At this point, we observe that the resulting effective 
Hamiltonian is in fact a Hamiltonian for the t-J model 
written in terms of the Jordan- Wigner fermions with the 
above constraint. By introducing the electron opera- 
acting in the restricted 
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Pit = a ]> Pol = aj/¥" Qj 
Hilbcrt space without double occupancies we obtain 



0(fc,w) = 
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Pkf\i>), 



(16) 



with the t-J Hamiltonian 



H=-t {phPja+h.c.)+J ^{S t ■ Sj+jmrij), 



(i,i> 



(17) 



where rij = nj pc and the hopping parameter t is 
defined as t = J/2. The ground state is now the 
tensor product of a vacuum state for holes and the \i/js) 
state. We have thus mapped the problem of a single or- 
bital excitation in the 2D FO and AF ground state, with 
dynamics governed by Hamiltonian Eq. (1) in the main 
text, onto a single hole doped into the undoped 2D AF 
ground state with its dynamics governed by Hamiltonian 

(LTZD- 



